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A computational approach is used to describe the aerothermodynamics of the Mars Pathfinder vehicle entering 
the Mars atmosphere at the maximum heating and maximum deceleration points in its trajectory. Ablating and 
nonablating boundary conditions are developed which produce maximum recombination of C0 2 on the surface. 
For the maximum heating trajectory point, an axisymmetric, nonablating calculation predicts a stagnation-point 
value for the convective heating of 115 W/cm 2 . Radiative heating estimates predict an additional 5-12 W/cm 2 at 
the stagnation point. Peak convective heating on the afterbody occurs on the vehicle’s flat stem with a value of 
5.9% of the stagnation value. The forebody flow exhibits chemical nonequilibrium behavior, and the flow is frozen 
in the near wake. Including ablation injection on the forebody lowers the stagnation-point convective heating 18% . 


Nomenclature 

Cj = mass fraction of species i 

D, = effective diffusion coefficient for species /, rrr/s 

I, = molar rate at which species i impacts surface, 

kg- mole/m 2 • s 

J, = mass flux of species i due to diffusion, kg/m 2 • s 
k r = surface reaction rate for reaction r, m/s 

Mt - species molecular weight, kg/kg-mole 
m a = ablation blowing rate, kg/m 2 • s 
P = pressure, N/m 2 

q w = convective heating at wall, W/cm 2 

R = universal gas constant, 8314.3 J/kg-mole K 
Ruxti ~ vehicle's effective nose radius, m 
s = distance along body surface, m 

T w = wall temperature, K 

V x = freestream velocity, m/s 

x, z = spatial coordinates, m 

v, = mole fraction of species i 

y, = participation fraction of gaseous reactant 

e - surface emissivity 

r) - spatial coordinate normal to wall, m 

p, = density of species /, kg/m 3 

<7 ~ Stefan-Boltzmann constant 

Introduction 

T HE Mars Pathfinder Project 1 - 2 proposes to land an unmanned 
probe on Mars to observe the planet’s surface and atmosphere. 
The project focuses on entry, descent, and landing as critical tech- 
nical challenges. To design the thermal protection system (TPS) for 
the entry vehicle, an accurate description of the aerothermal envi- 
ronment about the forebody and afterbody is necessary. 

Several recent studies 110 have predicted the,, aerothermody- 
namics on forebodies of vehicles entering the Mars atmosphere. 
References 3-5 address the Pathfinder mission directly. If the sci- 
ence payload is to be efficiently protected, however, the aerother- 
modynamics of the afterbody, or near-wake region, must also be 
described. 
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Computational methods can be used for this task. The thermo- 
chemical nonequilibrium Navier-Stokes tool described in Ref. 3 is 
used in the present work. It is an application of the LAURA (Langley 
aerothermodynamic upwind relaxation algorithm) of Gnoffo 11 12 to 
Mars atmosphere (assumed 97%' C0 2 and 3%; N 2 ). For the present 
study, the code has been further modified to perform an energy 
balance at the surface and include a fully catalytic C0 2 recombi- 
nation wall boundary condition with and without ablation-product 
injection. 

Wake fiowfields of acroshell-capped vehicles in air have been 
examined by numerous authors. Gnoffo has performed such a study 
with LAURA on lunar return and Mars return vehicles entering 
Earth’s atmosphere. 13 In addition, LAURA solutions have compared 
well with heating levels measured on the Shuttle Orbiter ieeside. 14 

The objective of the present study is to examine the aerother- 
modynamics about the vehicle entering the Mars atmosphere at the 
maximum heating point and the maximum deceleration point in its 
trajectory. Using results from a material response code for ablator 
response, a full (lowfield solution with ablating boundary condi- 
tions is also computed for the maximum heating point. Particular 
attention is given to estimating afterbody heating and describing the 
near-wake flow. 

Flowfield Computational Tool 

LAURA is an upwind-biased, point-implicit relaxation algo- 
rithm 12 for obtaining the numerical solution to the governing 
equations for three-dimensional, viscous, hypersonic flows in ther- 
mochemical nonequilibrium. The algorithm is based on a finite- 
volume formulation in which the inviscid fluxes are described by 
Roe’s averaging w ith second-order corrections based on Yec's sym- 
metric total-variation-diminishing scheme. The governing equations 
include continuity equations for 14 species, three momentum equa- 
tions (one of which is extraneous for the axisymmetric cases), and 
two energy equations describing the conservation of vibrational- 
electronic and total energies. For the present Mars atmosphere study, 
the 14 species are C0 2 , CO, N 2 , 0 2 , NO, C, N, O, SiO, H 2 , C 2 . Cj. 
H, and Si. (The last six species maintain trace concentrations ex- 
cept in the solution with ablation.) A 24-rcaction chemical-kinetics 
model is used. The reactions and rates are tabulated in Ref, 15. No 
appreciable degree of ionization is expected at these conditions. All 
transport properties for the species come from collision integrals, 
which are computed from Lennard-Jones potentials. Description of 
the transport properties and details on the remainder of the physical 
models are included in Refs. 3 and 1 1. 

Computational simulations of Mars entries suffer many limita- 
tions. The flow field is assumed to be steady and laminar. Unsteady 
effects are small at these high speeds, but transition to turbulence 
could occur, particularly at the maximum-deceleration trajectory 
point. Unfortunately, no available turbulence models have been 
shown to be valid in the wake regions of blunt bodies. In addition, 
many of the physical models used in the calculation are derived from 
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experimental data from compressive flows. Their applicability in the 
expanding regions of these flows is uncertain. Two additional limita- 
tions of the present method can produce elemental separation in the 
boundary layer. First, the multicomponent averaging procedure used 
to define the diffusion coefficients does not force the mass flux due 
to diffusion to sum to zero. Second, implementation of the flux lim- 
iters in the boundary layer may result in different difference stencils 
being used on different equations. Finally, transport properties may 
not be accurate at temperatures above 10,000 K or below 1000 K. 

Results presented are computed on a Cray C-90. An axisymmetric 
solution requires 3 to 6 h of computational time. 

Surface Boundary Conditions 

Species Continuity Equations 

In the present work, both ablating and nonablating boundary con- 
ditions are examined. In each case, since the degree of wall catalysis 
in C0 2 is unknown, the wall’s effect on the flow is assumed to cause 
maximum recombination of C0 2 Maximum recombination of C0 2 
should produce an upper-bound prediction for the convective heat- 
ing, and for this study, no other wall recombination is considered. 

Recombination of a dissociated gas mixture on a solid involves 
numerous complex phenomena. The present approximation esti- 
mates the upper limit of recombination by recognizing that such a 
process is limited by three rates: the rate by which reactants dif- 
fuse toward the wall, the frequency at which the reactants strike the 
wall, and the rate at which products can diffuse from the wall. Com- 
ment on the boundary condition’s extension to the ablating case is 
included at the end of the section. 

Assuming that the Rideal-Eley mechanism 16 is the only reaction 
mechanism, C0 2 surface recombination can be written as four first- 
order reactions: 


o + M -+ o , 

(1) 

CO + Oj — ► C0 2 + L] 

(2) 

CO + [, ] -»> CO, 

(3) 

0 + co., -*• co 2 + M 

(4) 


where the subscript .? refers to the surface-adsorbed reactant and [,] 
is an available adsorption site on the surface. 

The mass flux due to the reactions, J h must be balanced by diffu- 
sion. Using the same quasi-binary-diffusion approximation to multi- 
component diffusion used in LAURA, 11 we have 


dyo 

Jo = -^iPo - k 4 p 0 = ~pD 0 — 

dr) 


(5) 


Jco = ^2 Pco ~ hPco — -pDc o 


dyco 

drj 


( 6 ) 
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This approximation to multicomponent diffusion can be formulated 
in terms of the mass-fraction or mole-fraction gradient. The def- 
inition of the effective binary-diffusion coefficients is dependent 
on this choice. In the present work, as in LAURA, mole-fraction 
gradients are used, following the method of Lee. 1 

The reactions rates are 



/ RT W 

(8) 

1 In Mo 

*2 = K2y 

1 RT W 

(9) 

2n Mq 0 

11 

j RT W 

(10) 

2n Mco 

*4 = 

1 RT W 

(11) 

/ 2nM 0 


where y r is the fraction of gaseous-reactant particles striking the 
wall that participate in the reaction. For example, is the fraction of 
O atoms striking the surface that are adsorbed. To obtain a boundary 
condition from Eqs. (5-11) requires additional equations relating 
the Yr- 

Conservation of mass in a steady flow requires the surface cov- 
erage of O, and CO, to remain constant; thus. 



where c J UJ is the mass fraction of species i at the wall. With these 
relations, Jo+Jco+Jco 2 = 0 in Eqs. (5-7). Equations (12) and (13), 
along with the restriction that 0 < y r < 1, y\ + Yi 5 1. and 
y 2 + y 3 < L establish the boundary condition. In particular, if 


Cq.uj 



(14) 


maximum C0 2 recombination is caused by setting 


and 


Conversely, if 


yi + y4 

= 1 

(15) 

c 0.u. 1 

Y2 = Y\ 

CCO. 11 > 

( A/co\^ 
{mo) 

(16) 

Co, uj 1 
Y3 = Y4 j 

CcO, uj 

( A/co\* 
\M 0 ) 

(17) 

CcO.uj ^ Co, i/j | 

f Mco \ 5 
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(18) 


maximum C0 2 recombination is caused by setting 

y2 + y3 = i 


(19) 


and using Eqs. (12) and (13) for y ( and y 4 . 

A physical interpretation of these two cases can be extracted. 
From kinetic theory, the molar rate at which a gaseous species im- 
pacts the surface, /, , is the product of its local concentration with 
its thermal velocity divided by 4. For O, the expression is 


Po,w I 
~Mo V 2ttM 0 


For CO, the expression is 


PCO.u; / 

A/co V A/co 


( 20 ) 


( 21 ) 


If /o < A:o> then the inequality of (14) holds and the maximum re- 
combination of C0 2 is limited by the availability of O. In particular, 
the recombination in the reaction (2) cannot proceed faster than the 
maximum rate at which the reaction (1) can supply adsorbed O. In 
addition, since the CO, surface coverage is constant, CO cannot be 
adsorbed in the reaction (3) faster than the reaction (4) removes CO, . 

For the second case, when / C o < /o, the inequality (18) holds and 
the maximum recombination of C0 2 is limited by the availability 
of CO. 
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Expressions for species mass fractions at the wail can be written 
from Eqs. (5-7): 


held fixed in the boundary conditions, and the full vehicle flowfield 
is then computed. 


CO,w 


CO,w + 1 M m.w + 1 

M m<w + (k\ + ^ 4 ) (A 77 Mq/ D o) 


(22) 
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(24) 


The subscript w -f 1 refers to the location A r) above the wall. The 
wall temperature required in evaluating k { — k A is the radiative equi- 
librium value discussed below. The molecular weight of the mixture 
at the wall is obtained from the definition of molecular weight: 


Af m, w 


1 


(25) 


Geometry and Computational Grid 

A Viking-style geometry is proposed for the Mars Pathfinder entry 
vehicle. The forebody is a 70-deg sphere cone with nose radius of 
0.6625 m and an overall vehicle diameter of 2.65 m. A 0.06625-m- 
radius shoulder connects the forebody to a 49.7-deg cone afterbody. 
(Note: A longer vehicle is now being used for the Mars Pathfinder, 
whose afterbody cone is 47.44 deg. This slightly longer geometry 
is not examined in the present paper.) 

Creating a grid about the entire entry vehicle is a challenging 
problem. The hyperbolic grid generator HYPGEN 20 with HGUI 
(HYPGEN Graphical User Interface) is used to create the single- 
block grid. 

To assure accurate heating predictions, the grid is everywhere 
orthogonal to the body at its surface. The height of the first cell off 
the wall is 1.0 x 10“ 5 m, which yields a cell Reynolds number of 
3 based on the local speed of sound at the stagnation point. The 
grid extends 7.5 vehicle radii downstream. Surface grid distributions 
are generated using the GRIDGEN 21 22 software. Outer boundary 
limits are computed from previous elliptically generated grids. The 
orthogonality is controlled by HYPGEN in combination with three 
dissipative parameters. 


where the mass fractions for the remainder of the species are ob- 
tained from a noncatalytic wall boundary condition 

C/.u> — C|,u» + i ( 26 ) 

To preserve the simple explicit formulation for the mass fractions of 
Eqs. (22-24) in the numerical implementation, the value for the wall- 
mixture molecular weight is computed from the previous iteration’s 
wall mass fractions. 

Combining Eqs. (15-17) with Eqs. (22-24) yields equations for 
the CO, O, and C0 2 mass fractions at the wall when / 0 < lc o* Al- 
ternatively, when lc o < /o, Eqs. ( 1 2), ( 1 3), and ( 1 9) combined with 
Eqs. (22-24) provide the expressions for the wall mass fractions. 

When ablation injection is included, the blowing rate m a and 
pyrolysis gas composition c, 4 from a material response analysis 
appear in Eqs. (5-7) as 

dv; 

pD, — = m a (c i w - c Ls ) - J i (27) 

dr] 

Energy Equations 

Boundary conditions for the two energy equations require spec- 
ification of translational and vibrational-electronic temperatures at 
the wall. If the assumption of local thermal equilibrium is made at 
the wall, the two temperatures are equal. For the nonablating bound- 
ary condition, an energy balance at the surface sets the temperatures 
to the radiative equilibrium wall temperature. If q w is the convective 
heating at a location on the surface, then the wall temperature T w is 
obtained from 

q* = «7; 4 (28) 

where a is the Stefan-Boltzmann constant. The surface emissivity 
€ of SLA-561 char is assumed to be 0.78 on the forebody and 0.58 
on the afterbody. For the ablating boundary condition, the forebody 
wall temperature is fixed at the value predicted by the material- 
response analysis. 

Ablation Solution 

One solution is included with ablation injection defined by a 
coupled material-response analysis of the forebody. 18 The flowfield 
code employed in that coupled study used the same chemical kinet- 
ics, transport properties, and boundary conditions as those used in 
the present work. The material-response code was the CM A code 19 
working on the ablator material SLA-561 . In the present study, the 
blowing rate, composition, and wall temperature from the material- 
response analysis at the maximum-heating point in the trajectory are 


Results 

Based on a ballistic coefficient of 55 kg/m 2 and a 7.65-km/s ballis- 
tic entry at — 14.2 deg, one possible pathfinder trajectory is predicted 
to encounter maximum convective heating at 40.7-km altitude and 
6592-m/s velocity. The density and temperature of the C0 2 -N 2 gas 
at this altitude are 3.24 x 10~ 4 kg/m 3 and 162 K. This trajectory 
point is examined first and referred to as the maximum-heating 
point. The maximum-deceleration or maximum-pressure point ex- 
amined occurs at 28.5-km altitude at a speed of 4862 m/s. The 
atmospheric density and temperature are 1.161 x 1(T 3 kg/m 3 and 
176.8 K at this altitude. Flight Reynolds numbers based on vehicle 
diameter for these two trajectory points are0.64x 10 6 and 1 .53x 10 6 , 
respectively. 

Maximum-Heating-Point Results 

The forebody pressure prediction from a nonablating solution 
at the maximum heating point is compared with that predicted by 
a viscous shock layer (VSL) method 23 in Fig. 1. Both solutions 
display blunted cone distributions. The variable s is the distance 
along the surface from the stagnation point. An 5 value of 1 .437 m 
corresponds to the vehicle’s shoulder. Differences in the shoulder 
region might originate from differences in the two methods or from 
the VSL method’s forebody-only computation. 

The forebody convective heating predictions (nonablating) are 
compared in Fig. 2. The stagnation-point value from the present 
method is 1 15 W/cm 2 . whereas the VSL method predicts a value 
of 126.5 W/cm 2 . Both solutions predict a stagnation-point value for 
the C0 2 mass fraction of 0.83. In addition to different governing 
equations and solution techniques, the two methods use different 
diffusion models. 



Fig. 1 Comparison of pressure distributions on the forebody surface 
(V = 6592 m/s, alt. = 40.7 km). 
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Fig. 2 Comparison of convective-heating distributions on the nonab- 
lating forebody surface (V = 6592 m/s, alt. = 40.7 km). 



Fig. 3 Pressure and convective-heating distributions on the afterbody 
surface (V = 6592 m/s, alt. = 40.7 km). 



Fig. 4 Comparison of heating predictions on forebody and afterbody 
from grid resolution study in Ref. 15 (nonablating, V = 6500 m/s, alt. = 
40.5 km). 



Fig, 5 Close-up of streamlines in near-wake region (nonablating, V = 
6592 m/s, alt. = 40.7 km). 


Numerous approximate methods exist to estimate the stagnation 
heating on entry vehicles. The cold-wall correlation due to Sutton 
and Graves 24 for the Mars atmosphere is 

q u , = 1.89 x 10- ,,S /?-° e ^ 5 V^ (29) 

where is the effective nose radius in meters and is in W/cm 2 . 
When the actual nose radius of 0.6625 m is used, the correlation 
predicts a value of 1 19.7 W/cm 2 for this case. 

On the afterbody (Fig. 3), the highest pressure and convective 
heating occur on the flat stern. An s value of 2.83 m corresponds to 
the corner between the afterbody cone and flat stem of the vehicle. 
The maximum base pressure is 154 N/m 2 , or 1.1% of the stagna- 
tion pressure. The maximum base heating is 6.7 W/cm 2 , or 5.9% 
of the stagnation value. The vehicle’s drag coefficient is predicted 
to be 1 .69. 

Flight data from the two Viking missions 25,26 indicate that the 
base pressure is approximately 1% of the stagnation-point pressure 
for similar points in its trajectory. Preflight Viking predictions an- 
ticipated a base heating rate equal to 2% of the stagnation-point 
value. Reference 25 states the flight data exceeded that prediction 
and reports a value of 4.2% for apparent base heating. The Viking 
probes, however, entered at 1 1 -deg angle of attack and speeds near 
3500 m/s. In addition, the locations of the pressure transducers and 
thermocouples were not necessarily the locations of maximum pres- 
sure and heating. 

A grid resolution study was performed in Ref. 15. The freestream 
velocity of 6500 m/s and density of 3.19 x 10~ 4 kg/m 3 for that 
study are similar to the maximum-heating conditions examined here. 
Figure 4 (from Ref. 1 5) compares the convective heating prediction 
from a 136 x 80 grid like the one used in the present study with 
a 136 x 160 grid. Based on this comparison, the 136 x 80 grid is 
judged as adequate. 

In the present axisymmetric solution, the maximum afterbody 
base pressure and heating are the result of impingement by a strong 
recirculation vortex. A close-up of this vortex is shown in Fig. 5. The 



Fig. 6 Contours of CO 2 mass fraction (nonablating, V = 6592 m/s, alt. 
= 40.7 km). 

vortex extends downstream to a free- shear- layer stagnation point 
5.17 m from the nose. A weak secondary vortex exists just behind 
the shoulder and is also indicated in Fig. 3. (The three-dimensional 
manifestations of these axisymmetric vortices are toroidal in shape.) 

The flow in the nose and conical regions of the forebody is 
subsonic. The subsonic-zone deflection angle around the vehicle’s 
shoulder is 1 1 deg. 

A question arises as to the chemical state of the near wake. Rapid 
expansion of reacting flows can result in chemically frozen flow. 
Figure 6 presents contours of C0 2 mass fraction for the near-wake 
flowfield. The figure indicates the source of C0 2 in the wake results 
from the boundary condition, which forces C0 2 recombination. The 
C0 2 so recombined on the forebody and afterbody is then entrained 
in the recirculation vortex. 

The radiative heating expected for the maximum-heating case is 
estimated by several methods. The equilibrium method of Ref. 27 
predicts a value of 5.5 W/cm 2 at the stagnation-point. A more de- 
tailed forebody calculation using the LORAN code of Hartung 28 
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to that shown in Figs. 5 and 6, except the primary recirculation vor- 
tex is smaller, with the free-shear-layer stagnation point at 4.41 m 
downstream of the nose. The drag coefficient for this case is pre- 
dicted to be 1.71. The stagnation-point radiative-heating estimate 
for this case is less than 0.1 W/cm 2 . 

Conclusions 

The convective heating at the stagnation point computed for the 
nonablating maximum-heating trajectory point is 1 15 W/cm 2 . The 
maximum afterbody heating and pressure for this zero-angle-of- 
attack case occur on the vehicle’s stern. These maxima result from 
the focused impingement of a large recirculation vortex, and the 
maximum afterbody pressure is near 1% of the stagnation pressure. 
The maximum afterbody convective heating is 5.9% of the stagna- 
tion value. These ratios are consistent with Viking flight data. Ra- 
diative heating on the forebody is 5-11% of the convective heating. 

The effect of forebody ablation with blowing rate at 1% of 
freestream mass flux is to decrease the forebody convective heating 
18%' at the stagnation-point. The decrease on the forebody flank is 
as high as 50%. Products from the forebody ablation become en- 
trained in the recirculation vortex, but they have little effect on wake 
structure and afterbody surface pressure and heating. 

Nonablating results from the maximum deceleration point in 
the trajectory estimate the maximum stagnation-point pressure at 
26,525 N/m 2 . While forebody convective heating levels are 40-50% 
lower than those predicted at the maximum-heating trajectory point, 
the afterbody convective-heating values are nearly equal. 
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